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Abstract 

We propose to study factorization breaking effects in exclusive b decays where they are 
strongly enhanced over the factorizing contributions. This can be done by selecting 
final-state mesons with a small decay constant or with spin greater than one. We 
find a variety of decay modes which could help understand the dynamical origin of 
factorization and the mechanisms responsible for its breaking. 
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1 Introduction 



An outstanding task in heavy-flavor physics is to understand the strong interaction effects in 
exclusive weak decays of hadrons containing a 6-quark. For many decay channels, such an 
understanding is a precondition for gaining information on the quark mixing matrix or on 
physics beyond the standard model. In addition, the dynamics of quarks, gluons, and hadrons 
in the presence of a large mass m& is interesting in QCD by its own right. 

Introduced in U, the concept of factorization has been one of the most successful tools in 
this respect, providing fair agreement between theory and data for many channels. In other 
cases, factorization in its most naive version fails when compared with experiment, and there 
have been several phenomenologically motivated improvements over its original form Q. || . 

There are several dynamical arguments why and where factorization should be valid. One 
is based on the large N c limit of QCD whereas a different line of approach builds on 
the color transparency phenomenon 0]. More recently, the framework of QCD factorization 
has implemented the color transparency argument in the language of perturbation theory and 
power counting in l/m b || |7|, ||. 

In these approaches it is understood that there are corrections to naive factorization, which 
are suppressed in a small parameter such as 1/N C , or a s and AQ CD /m fe . A more quantitative 
understanding of their size is crucial in order to assess for which channels and to which precision 
the factorization concept can be applied. There are also scenarios where factorization in the 
sense of does not appear as a limit when a small parameter vanishes, and where conceptually 
factorization breaking terms in the decay amplitude may be as large as the factorizing ones. 
An example is the perturbative hard scattering (PQCD) approach ||. In view of such contro- 
versies, and given that present day theory can at best estimate the size of most nonfactorizing 
contributions, quantitative tests of factorization in the data are of great importance. 

We propose here to study decay channels where the factorizing contributions to the ampli- 
tude are small or zero for symmetry reasons. In such a situation nonfactorizing contributions, 
which would otherwise be suppressed, have a chance to be clearly visible. The measurement 
of the corresponding decay rates can thus give rather direct information on their size, and the 
comparison of different channels may give indications on the relevant dynamical mechanisms. 

Our suggestion is to choose decay channels whose flavor structure is such that a selected 
meson X must be emitted from the weak current mediating the 6-quark decay. Taking then a 
meson which has as very small decay constant, the factorizable contributions to the decay are 
suppressed. A second possibility is to consider mesons X with spin J > 2. A tensor meson 
for instance cannot be produced from a decaying W boson, which has spin 1, unless there 
are interactions involving the other hadrons in the decay process. We will find a variety of 
decay channels where these ideas can be realized, which will allow us to address different issues 
related to factorization and its breaking. 

The organization of this paper is as follows. In Sect. we review the basics of factorization 
which will be essential for our arguments. We select the mesons for which factorizing contri- 
butions in decays are suppressed in Sect. |3|. In the following section we identify which flavor 
structure a decay must have in order for this suppression to apply, and take a closer look at 
specific issues in the different channels. Sect. |5| discusses how suppression can be circumvented 
by different nonfactorizing mechanisms. Some of these can be treated within QCD factoriza- 
tion and will be investigated in Sect. We estimate branching ratios of suppressed decays 
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into a heavy and a light meson in the subsequent section, before concluding in Sect. Some 
numerical estimates concerning meson distribution amplitudes, which we need in our paper, 
are given in an Appendix. 



2 Matrix elements for hadronic two-body decays 



We start by briefly recalling the low-energy effective Hamiltonian and some basics of the factor- 
ization approach. Hadronic two-body 6-decays are described by the effective weak Hamiltonian 



nj G F 



£ v jb v k * d (c^ + c 2 oi k ) - v tb v td £ ca 

j,k=u,c i 



{d — > s} + h.c. 



1) 



where V denotes the CKM matrix. The operators 0{ k 2 result from tree level W exchange and 
in the case i = k = u read 



QUU 

W 2 



w a 7 M (l - 7 5 )6 Q dpi^l - 75)11/9, 



(2) 



Here a and j3 are color indices, and it is understood that all fields are taken at space-time ar- 
gument zero. The remaining operators in H e s are so-called penguins. For a detailed discussion 
of the operators Oi and the Wilson coefficients Cj we refer to 
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In naive factorization, the matrix element (YX\ H. c g \ B) is written as a product of matrix 
elements of quark currents between B and Y, and between the vacuum and X. Only the 
color-singlet piece of each current is retained, while the color-octet piece is neglected. This 
leads to replacing 7i cS by the effective transition operator T, whose tree level operators read 



T {1 > 2) = ^V ub V: d [a 1 uni-K)b®d 1 ,(l- l5 )u 
+ a 2 d^(l - 75)6 ® «7 M (1 - 7 5 )w 



(3) 



for i = k = u. Here the notation ® indicates that the matrix elements are to be taken in 
factorized form as described above. The new coefficients a± = C\ + C 2 /3 and a 2 = C 2 + Ci/3 
have been obtained from projecting on color-singlet currents, and are commonly referred to 
as color allowed and color suppressed, respectively. Numerically, a\ is close to 1 and a 2 of 
order 0.1 at a renormalization scale \i = m b . Notice that the Fierz transform performed in the 
second term of Eq. (^) has left the (V — A) x (V — A) structure invariant, where V and A 
respectively denote the vector and axial vector current. The situation is analogous for those 
penguin operators which again involve (V — A) x (V — A) currents, for explicit formulae see 
e.g. @. 

Important for us will be the strong penguins with (V — A) x (V + A) structure (the 
electroweak ones with similar Dirac structure are numerically less important in the Standard 
Model). They are 



5 = d a Y(l ~ lb)b a It^ni 1 + 75) 9/?, 

6 = 47 M (1 - 75)&/3 E<?/37m(1 + 75)<? a , 



(4) 
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X a b x 7r a 2 a tt 2 p 3 vr a 4 Xco Xc2 

m x [MeV] 985 1230 1300 1318 1474 1670 1691 1801 2014 3415 3556 
jpc ++ 1+ - -+ 2 ++ ++ 2 -+ 3— 0-+ 4++ 0++ 2++ 



X 


K* 


K* 


^2 


K 3 * 


^2 




^2* 


P s j 


m x [MeV] 


1412 


1426 


1773 


1776 


1816 


2045 


2459 


2573 


J p 


0+ 


2+ 


2~ 


3~ 


2" 


4+ 


2+ 


(2 + ) 



Table 1: Mesons for which one of the suppression mechanisms discussed in Sect. is relevant. 
Masses are taken from the Review of Particle Data 
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We do not list those mesons that are 
omitted in the meson summary table. The spin-parity assignment for the D s j(2537) is not 
certain yet. 



with a sum over q = u,d, s, c, b, and in naive factorization lead to 

G p 



7-(5,6) 



^2 VtbVtd 



+ a 6 J2 (- 2 ) 9(1 - 7s)& ® d(l + l 5 )q 



(5) 



The structure (P — S) x (P + S) involving the scalar and pseudoscalar currents S and P has 
emerged from the Fierz transform in the term with = Cq + C^/3, whose value is about —0.03 
at \i — m&. The corresponding operator provides one possibility to circumvent the suppression 
mechanisms we will discuss shortly, and we will often refer to it as scalar penguin. 



3 Meson candidate selection 

Let us now specify our mechanisms to suppress decay amplitudes in naive factorization and 
see to which final states they apply. 

3.1 Suppression mechanisms 

There are several reasons why the coupling of a meson to the local currents of the effective 
weak Hamiltonian can be suppressed. Clearly, a meson with spin J = 2 or larger has no matrix 
element with either of the currents S, P, V, A, and in naive factorization cannot be produced 
as a meson ejected by the effective weak current. From Table [1] we see that examples for such 
mesons are the a 2 , vr 2 , p 3 , and K\ in the light quark sector. Heavy tensor mesons are the D\ 
and D s j, and there is a tensor charmonium state, \2c- 

Let us now turn to mesons with J = 0, 1 whose production from the weak current in naive 
factorization is forbidden or suppressed because their coupling to V and A is zero or small. 
We define the decay constants of the negatively charged mesons with isospin I = 1 as 

(S(q)\d(0)Yu(0)\0) = -if git, (6) 
(P(q)\d(0) 7 M 75«(0)|0) = -i/pg", (7) 
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(V(q,e)\ J(0)7 M «(0)|0) = -if v m v f, 
(A(q,e))\ J(0) 7 M 75«(0)|0) = -if A m A e», 



for scalar, pseudoscalar, vector, and axial mesons, respectively. Here denotes the meson 
momentum and, if applicable, & its polarization vector. The choice of vector or axial quark 
current on the right-hand sides of these definitions is dictated by the parity of the meson. For 
the corresponding neutral mesons, the flavor structure of the current is (uu — dd) / v2 instead 
of du, and for mesons with different quark content one has to take sd, su, cc, etc. 

Because of charge conjugation invariance, the decay constants for the neutral ao and b\ 
mesons and for the Xoc must be zero. In the isospin limit, the decay constants of the charged 
a and b\ must thus vanish, too, so that f a0 and foi defined in Eqs. (|), @ are small, of order 
m d — m u- For the a mesons, this can explicitly be seen by taking the divergence of Eq. (|]), 
which by virtue of the equations of motion gives 

<J ao = i(m d -m u )(a o \d(0)u(0)\0). (10) 
For the charged Kq the analogous relation reads 

m 2 K *f K * = i(m s -m u )(K*\s(0)u(0)\0), (11) 

which becomes zero in the flavor SU(3) limit and indicates that f%* should be suppressed. 
It is instructive to compare Eqs. ( |10|) and (|Tl|) with their analogs for the pseudoscalars, 

m 2 J n = i(m d + m u ) (tt| J(0) 75 «(0) |0), (12) 
m 2 K f K = i(m s + m u )(K\s(0) l5 u(0)\0). (13) 

Because the axial current is not conserved, these decay constants do not vanish in the isospin 
or SU(3) limit. Moreover, they do not vanish in the chiral limit, m u = = m s = 0, for 
the light pion and kaon, since these mesons are Goldstone bosons and become massless in the 
same limit. Numerically, f w and fx are in fact not small and of the same order of magnitude 
as for instance f p and fx*- The decay constant for the heavy 7r(1300), however, does become 
zero in the chiral limit, and its actual value is small due to chiral suppression. 

We remark that the spin-zero mesons whose coupling to the V and A currents is small or 
zero for one of the above reasons can still couple to the S or P currents appearing in penguin 
operators of the effective Hamiltonian, as discussed in Sect. |2|. This does however not hold for 
the b\, which has no matrix elements with S or P. 

3.2 Decay constants 

The decay constants of the a (980), a (1450), Kq, 7r(1300), and the b\ are poorly known at 
present. Using finite energy sum rules, Maltman 113] obtained^ 

/ a0(980 ) = 1.1 MeV, / a0{1450) = 0.7 MeV, f K * = 42 MeV, (14) 

1 Maltman defines the ao decay constants with an extra factor (m s — m u )/(md — m u ). To convert them into 
our convention, we take the quark masses in Eq. dig) below. 
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consistent with the ranges estimated by Narison [T3j 

iaooso) = 0.7 to 2.5 MeV, f K . = 33 to 46 MeV. (15) 
For the heavy pion, the theoretical estimates in [|14j] provide a range 

A(i300) = 0.5 to 7.2 MeV. (16) 

Comparing these values to 

f n = 131 MeV, f K = 160 MeV, (17) 

we find that the suppression patterns discussed in the previous subsection are indeed seen 
numerically, with the decay constants for the ao mesons smaller than those for 7r(1300) because 
of the relative signs between the quark masses in Eqs. (|TI]) and (|TJ). We also see that f K * 
is suppressed relative to but not as strongly as f ao compared with f n , because SU(3) 
symmetry breaking is rather strong for the quark masses. Here we have implicitly assumed 
that the (pseudo) scalar matrix elements on the right-hand sides of Eqs. fllOD to (|1^) are not 
anomalously small or large. Indeed, taking from fl5| 

m u = 4.8 MeV, m d = 8.7 MeV, m s = 164 MeV (18) 

for the MS quark masses at \i = 1 GeV, we find with the values (|16|) and ( |17|) 

i(rr\d-f 5 u\0) « 0.19 GeV 2 , 

i (tt(1300)| J75M |0) « 0.06 to 0.90 GeV 2 , 

i (K\ |0) « 0.23 GeV 2 , (19) 

and with <M) and (M 



z(a (980)| du\0) w 0.17 to 0.62 GeV 2 , 
i(ao(1450)|d«|0) w 0.39 GeV 2 , 

i (^(1430)1 su |0) w 0.42 to 0.58 GeV 2 . (20) 

Despite a certain spread these values are remarkably close to each other, given that the cor- 
responding squared meson masses vary by more than two orders of magnitude. We note that 
Chernyak |16|| has recently estimated fx* — (70 ± 10) MeV. We consider this to be rather high 



as it is far away from the range (O) obtained in other studies. Also, the corresponding value 
for (Kq (1430)| su |0) would correspond to a quite strong SU(3) breaking for the scalar matrix 
elements. 

We wish to emphasize at this point that the decay constants for the a (980), a (1450), Kq, 
7r(1300), and the b\ can be measured very cleanly in r decays. In fact, from the bound on 
the branching ratio B(t — > 7r(1300)z/ T ) < 1 • 10~ 4 in we infer ^(1300) < 8.4 MeV, which 
is not far from the upper end of the theory estimates (|I6|) . The decay constants in Eq. flTil) 
correspond to branching ratios of 

B{t -> a (980) u T ) ~ 3.8 ■ 10" 6 , 
B{t -> a (1450) v T ) ~ 3.7 • 10" 7 , 

B(t -> K* (1430) v T ) ~ 7.7 -lO" 5 . (21) 



2 We have taken the average values given in Table 6 and evolved them down using Eq. (20) in that reference. 
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0.5 1 1.5 2 2.5 

m x [GeV] 



Figure 1: Values of E x /m x as a function of m x for the decays B — > DX (lower curve) and 
B — > 7iX (upper curve). 



These estimates are rather encouraging, given that at the B factories one expects to have 
about 3 • 10 7 t pairs with 30 fb ]17| , and that the potential of r-charm factories would be 



even higher. A measurement of the decay constants for the above mesons would greatly reduce 
the uncertainties in the predictions we will give in Sect. [7[ It could also provide valuable 
information on the nature in particular of the ao(980) and ao(1450), only one of which can be 
a member of the conventional qq meson nonett. 



3.3 Kinematics 

The color transparency argument || for factorization of decays B — > YX requires the meson X 
emitted from the weak current to be fast. More quantitatively, its time dilation factor Ex/mx 
should be large, where Ex = {m 2 B — m\ + m x )/(2rriB) is the energy of X in the B meson rest 
frame. We show the values of Ex/mx for Y = D and Y = it in Fig. |l|. The corresponding 
curves for B — > D*X and B s — > D^*'X are very close to the one for B — > DX, and the ones 
for B — > pX, B s — > KX, and B s — > K*(892)X are practically the same as for B — > ttX. 
Only if X is a pion does one have a very large Ex/mx, namely E n /m n = 16.5 for B — > Dtt 
and E^/m^ = 19 for B — > tttt. For mx above 1 GeV, relevant for the mesons in Table [I], this 
ratio decreases rather gently from moderate values down to little above 1. Even so, the mass 
range of our candidate mesons seems sufficiently large so that a study of the corresponding 
decay channels could provide valuable clues to whether corrections to factorization significantly 
depend on Ex/mx, or more generally, on the mass mx- We remark in passing that even the 
lowest value of Ex/mx in Fig. [l| corresponds to a velocity fix = 0.5 and a recoil momentum 
of px = 1.5 GeV in the B rest frame, indicating that X is still relativistic. 
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3.4 Resonance decays and the continuum 



Clearly, the measurement of rare B decays involving higher mass resonances presents an experi- 
mental challenge. An experimental analysis will be easier if the meson X in question has a decay 
channel with sizeable branching fraction that is not also accessible to mesons with nearby mass 
whose production is not suppressed. To give some examples, the decays of a (980), a 2 (1320), 
and ao(1450) into 7irj appear rather clean in this respect, as do the modes &i(1235) — > urn and 
tt 2 (1670) -> / 2 (1270)tt. On the other hand, the decays of a 2 (1320), vr(1300), and vr 2 (1670) into 
pir are more problematic because of background from the rather broad cti(1260). The same 
holds true for the decays of the _K"o(1430) and 7^2(1430) into Ktt because of background from 
the AT* (1410). In such partial wave analysis of the decay products will probably be 

necessary in order to constrain the decays into the mesons with spin J = or J = 2. 

We emphasize at this point that mesons X with a rather large decay width do not present 
as serious a problem in our context as in other studies. The physical arguments leading to 
a suppression of their production within the factorization mechanism (as the arguments for 
factorization itself) do in fact not depend on X being a narrow resonance. Our arguments in 
Sect. |Oj were at the level of current matrix elements and go through in complete analogy if, 
for instance, the |7r(1300)) state in Eq. (|T2| ) is replaced with \p7r) in appropriate partial waves 
that have definite quantum numbers J p = 0~. Moreover, the main idea here is to use the 
branching ratios of suppressed decays as quantitative estimates for the size of corrections to 
factorization, and to study their pattern by comparing different decay channels. An uncertainty 
on the branching ratio of B — > YX due to the line shape of X is therefore less severe as in 
channels which are allowed by factorization, and where factorization tests need branching ratios 
to a much higher precision. 

One could in fact also perform the studies we propose here not with particular meson 
resonances but with continuum states, similarly to a recent test of factorization by Ligeti et 
al. ||18|| . Our main reason to concentrate on resonances X here is that, by definition, their 
production should be enhanced with respect to continuum states with the same quantum 
numbers, which is important since we are looking for decays with small branching ratios from 
the start. 



4 Decay mode selection 

We are looking for exclusive decays B — > YX, where the meson X must be emitted from the 
weak decay vertex and cannot pick up the spectator quark. Only then will the factorizing 
contribution to the decay be suppressed for mesons X with small or vanishing decay constant 
or with spin J > 2. This puts requirements on the flavor structure of the decay, which we 
now discuss. One may avoid these requirements by studying decays where both final state 
mesons are taken from Table [I], such as B° — > bfb± or B° — > . F° r one °f the mesons the 
appropriate suppression mechanism will then always be at work. 

For definiteness we consider in the following the case where the B meson contains a b and 
not a b. One requirement now is that the flavors of the spectator antiquark and of the antiquark 
emitted from the b decay must not be the same, otherwise X can pick up either of them. We 
thus cannot use decays such as B~ — > D°ciq , whose flavor structure reads ub — > u(cud), where 
the brackets indicate the quarks originating from the electroweak vertex. 
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Figure 2: Tree level diagrams for (a) B° — > D + a and (b) B° — > D + a, 



A second requirement is due to imperfect knowledge of the initial state. Whereas the decay 
B° — * D + dQ with flavor structure db — > d(cud) satisfies our requirement, the same final state 
can be reached in the decay of a 5°, where the ag contains the spectator. The corresponding 
diagrams are shown in Fig. |^. This CP conjugated background can in principle be removed 
by flavor tagging, which puts of course stronger demands on the experiment. In the example 
just given, the amplitude from B° decay is suppressed by A 2 relative to the direct decay, where 
A ~ 0.22 is the Wolfenstein parameter in the CKM matrix. We will give a more quantitative 



estimate of this background in Sect. 7A . In other cases however, e.g. for B s — ► D+Kq~ versus 
B s — > D+Kq~, both signal and background amplitudes are of order A 3 . We will not consider 
such modes in the following, and list in Tables ^| and [3] the flavor structure of decays satisfying 
the following conditions: 

1. It is ensured that the meson selected from Table [1] cannot pick up the spectator antiquark 
in the B meson. 

2. A CP conjugated background that violates condition 1 either does not exist or is CKM 
suppressed. It turns out that the only case we retain that has such a background is the 
one mentioned above, listed in the first row of Table 0. 

Let us now consider the different decay categories separately. 



4.1 Decays with one or two heavy mesons 

Decays into heavy-light final states with open charm are the simplest from the point of view 
of their electroweak structure, since only the W exchange operators 0\ and O-i of the effec- 
tive Hamiltonian contribute. For color allowed decays B — > DX with X = ir,p, a 1; naive 
factorization is in rather good agreement with data || 

For color allowed decays where the heavy meson is the emission particle, the color trans- 
parency argument does not hold, but arguments based on the large N c limit do. Relevant 
meson candidates here are the D\ and the D s j. Comparing the size of nonfactorizing contribu- 
tions for the two types of channels with the suppression mechanisms discussed here might thus 
shed light on the question which type of mechanism is more relevant to ensure factorization. 
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example decay 


factorizing 


contribution 


annihilation 




quark level 


tree peng. 


topology tree peng. 


B° -> D + a 


c?6 — > d(cud) 


A 2 


a A 2 




— ► d(cus) 


A 3 




B s - D+ao 


sb — > s(cug?) 


A 2 





£° -> tt+ZTj dfe -> d(ucs) A 3 

B- -> 7r Z>2 _ ub -> w(wcd) A 4 a A 4 

— > "^-D^ - "~ > u{ucd) A 4 a, c A 4 

— > 7r°D~j — ► u{ucs) A 3 

— > — * u{ucs) A 3 a, c A 3 

B s -> iT+Dr s6 -> s(wcd) A 4 







db- 


-> J(ccs) 


A 2 


A 2 


a 




A 2 


5- - 




ub - 


-> u(ccd) 


A 3 


A 3 


a 


A 3 


A 3 








-> -u(ccs) 


A 2 


A 2 


a 


A 4 


A 2 






sb - 


-> s(ccd) 


A 3 


A 3 


a 




A 3 






db- 


-> d(uus) 


A 4 


A 2 


a 




A 2 


5- - 




ub - 


-> u(dds) 




A 2 


a 


A 4 


A 2 








-> u(ssd) 




A 3 


a 


A 3 


A 3 




- X+a 2 - 


sb~ 


■> s{uud) 


A 3 


A 3 


a 




A 3 



Table 2: Color allowed decay modes -B — > FX satisfying the criteria specified in the text. 
X is the meson emitted from the weak decay vertex. Any meson can be replaced by another 
with the same flavor structure and isospin, e.g, D by D*, 7r by p, and the meson X by an 
appropriate candidate from Table [3]. The second column gives the quark composition of the 
final state, with the quarks originating from the b decay enclosed in brackets. In the third 
and fourth columns we give the power of the Wolfenstein parameter A in the decay amplitude 
for tree level W exchange and penguin operators, respectively. Annihilation contributions are 
listed with their topology as shown in Fig. |5] and the power of the Wolfenstein parameter for 
tree level and penguin operators. 
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example decay factorizing contribution annihilation 

quark level tree peng. topology tree pen£ 



B°- 


-> n°D? 




-> d(duc) 


A 2 


a 


A 2 








-* d(duc) 


A 2 


a, c 


A 2 


B~ - 


-»■ K~D*° 


ub - 


-> m(smc) 


A 3 






B s - 


■+ fT /^ 


sb- 


-> s(duc) 


A 2 














A 3 


ct . C 


A 3 



B ^K Xc0 db^d(scc) A 2 A 2 b A 2 

^7r x c o ^d{dcc) A 3 A 3 b A 3 A 3 

^ VXco — ► rf(rfcc) A 3 A 3 b, c A 3 A 3 

B- -> i^x c o m6 -> w(scc) A 2 A 2 b A 4 A 2 

— ► 7r _ Xco — ► w(c?cc) A 3 A 3 b A 3 A 3 

B s ^ i] Xco sb — > s(scc) A 2 A 2 b, c A 2 A 2 

^K° Xc0 ^s{dcc) A 3 A 3 b A 3 



5 S - 




s6- 


-> s(duu) 


A 3 


A 3 


a 




A 3 








-> s(smt) 


A 4 


A 2 


£L ? C 


A 4 


A 2 



Table 3: As Table || but for decays where the tree level contribution is color suppressed. We 
remark that the decay B s — > K° a® has a color allowed penguin contribution. 

One may also use decays into two charmed mesons to address the same question. We remark 
in this context that in recent study, Luo and Rosner found factorization to work reasonably 
well for B° — > within present errors figfl . 

Notice that for color suppressed channels such as B° — > -k q D^ naive factorization is neither 
backed up by color transparency nor by 1/N C arguments. The decays into a in Table § 
will thus show whether the factorization concept still applies here. 

Let us finally consider B decays into charmonium. Naive factorization has notorious prob- 
lems with these channels ||, compounded by the fact that the coefficient 02 is extremely 
dependent on the factorization scale \i. A comparison of the decays involving a Xco 01 Xc2 
and the corresponding ones with J ftp or Xci m &y thus shed light on the relative importance of 
factorizable and nonfactorizable contributions. 

4.2 Penguins and decays into two light mesons 

Decays into two light mesons present some specifics due to the presence of penguin operators. 
First, we have no meson candidates with isospin 1 = 0, which are a superposition of quark 
states mm, dd, ss. Penguin transitions lead to all three of them, and one will always violate our 
condition that the spectator and the emitted antiquark must have different flavors. Second, as 
remarked at the end of Sect. |3.1| , suppression mechanisms based on the smallness of the decay 
constant fx are not effective if scalar penguins occur in the factorization ansatz. Only spin 
suppression, and the isospin suppression for the 61, are still at work then. 



11 



It is nevertheless instructive to look at the relative importance of the S and P oper- 
ators in decays involving scalars or pseudoscalars. In the decays B° — > ajj" dp or B° — > 
7r + (1300)7T _ (1300) the scalar penguins come with huge enhancement factors in the amplitudes 

r a = 2m l r -K = 2m l / 22 x 

Numerically, r n( - 130 ^ = 85, r a °( 98 °) = 170, and r a °( 145 °) = 380, to be compared with r n = 1 for 



the light pion, where we have evolved the light quark masses Q18|) up to /i = m& = 4.4 GeV. 
The strong scalar penguin can hence compete with the current-current operators, even though 
its coefficient introduced in Sect. ^| is only of order of several 10 -2 . 

Using Eqs. fllCf) , fll2|) , (p2]) , we can express the products fxT x in terms of the (pseudo) scalar 
matrix elements ( |i~9"D and (|20|), and observe that 

f ao r ao « A(i300)r^ (1300) « (23) 

i.e., they are all roughly of the same size. Since B° — > 7r + 7r~ is driven by the color allowed 
tree level coefficient a±, naive factorization predicts the decay rate for B° — > aQ to be small 
compared with the one for B° — > 7r + 7r _ . Namely, the ratio of their amplitudes is controlled by 
the small parameters f ao / f n and 

oe fa r ao = ^6 r n {ao\_du |0) ^ 24 ^ 

Ol /,r Oi (7r| J75 U |0) ' 

corresponding to the tree level and scalar penguin contribution to B° — > Oq a^, respectively. 
Analogous estimates can be given for B° — > 7r + (1300) 7r _ (1300), and also for B s decays into 
K + clq and i^ + 7r _ (1300) compared to K + n~. 

The situation is different for decays where the emitted meson is a kaon. The decay mode 
B° — > 7i + K~ is penguin dominated since its tree level contribution is CKM suppressed. With 
the analog of ( p3|) for strange mesons one thus obtains similar decay rates for B° — > 7r + K~ 
and B° — > tt + Kq~, where the latter receives most of its contribution from the scalar pen- 



guin, as was pointed out in [|16| . We emphasize that, in contrast, naive factorization predicts 

b(b° -> w+k;-) = 0. 

4.3 Bottom baryon decays 

Bottom baryons provide a complementary field to study exclusive hadronic decays, making 
more degrees of freedom such as polarization accessible to experimental investigation. The 
most notable differences between qq and qqq bound states in the context of factorization studies 
are the quark content and the role of annihilation topologies. Since the initial baryon can never 
be completely annihilated by the operators in Eq. (P, we call the corresponding topologies 
shown in Fig. ^ pseudo-annihilation. For an overview of heavy baryon decays, we refer to p0 |. 



Let us adapt our ideas to study factorization and its breaking with spin or decay constant 
suppression to the case of exclusive heavy baryon decays. Of course there is no background 
here from decays of the CP conjugated parent into the same final state, such as discussed in 
Sect. |4]. In order to ensure the formation of the final state meson from the electroweak current, 
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— > ud(ccs) 


A 2 


A 2 


A 4 


A 2 




— > ud(ccd) 


A 3 


A 3 


A 3 


A 3 



fit, — > fi c a , H -D2 ss b ~^ ss(cud) A 2 

-> VLDf ^ss(cus) A 3 

-> fii^ ^ss(ucs) A 3 

-> H ^-, -> ss(wcd) A 4 

— ► ^Xco — ► ss(ccs) A 2 A 2 A 2 

-> fi c ZAp> -~Xco — > ss(ccd) A 3 A 3 A 3 

-> Vtal -^ss{uus) A 4 A 2 

— ► H _ a° — > ss(uud) A 3 A 3 A 3 



Table 4: Flavor structure of the decay modes of the A& and the fib for which the emitted meson 
must originate from the weak current. The mesons can be replaced by others from Table [I] with 
the same flavor structure. The second column gives the quark composition of the final state, 
with the three quarks originating from the electroweak vertex enclosed in brackets. The third 
and fourth columns give the power of the Wolfenstein parameter A in the decay amplitude 
for W exchange and penguin operators, respectively, and the corresponding information for 
pseudo-annihilation contributions is given in the last two columns. 

we must however require that the spectator quarks in the baryon be different from the quarks 
produced in the weak decay. In addition, we can only consider bottom baryons that have weak 
decays and do not dominantly decay strongly or electromagnetically. A possible decay channel 
is for instance Ab —>■ A c D~j. Also, Mannel et al. have mentioned that the fib might only 
have electroweak decays. In Table [| we list the decays of these two baryons for which it is 
assured that the meson cannot pick up a spectator, so that corrections to factorization can be 
studied with our method. 



5 Escaping suppression by factorization breaking 

The idea of this paper is to study the size and pattern of corrections to factorization in an 
environment where they are not "hidden" behind a larger factorizing piece. Without giving 
an exhaustive discussion of nonfactorizing contributions, we now consider two of them and see 



why the various suppression mechanisms discussed in Sect. [3.1| do not apply. 
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Figure 3: The different topologies for pseudo- annihilation contributions to b baryon decays. 




(a) (b) 

Figure 4: Examples of diagrams with gluon exchange breaking naive factorization: (a) vertex 
correction, (b) spectator interaction. 

5.1 Nonfactorizing gluon exchange 

Naive factorization is broken by strong interactions of the quarks originating from the b decay. 
In the language of quarks and gluons they correspond to diagrams like the ones in Fig. |j. What 
is important in our context is that such contributions no longer involve the matrix elements of 
the meson X with the local currents V, A, S, P. As a consequence the suppression mechanisms 
based on the smallness of the decay constant fx are not effective. Also, one or several gluons 
absorbed by the quark-antiquark pair that will form X can transfer both helicity and orbital 
angular momentum, so that the spin of X is no longer restricted to be or 1. 

Note that these arguments are independent on whether the internal lines in the diagrams 
of Fig. [|have large virtualities or not. In the first case the corresponding contributions can be 
calculated in perturbation theory. In Sect. ^| we will analyze them in the QCD factorization 
approach and explicitly see that our suppression mechanisms are no longer operative. 

If the internal lines in these diagrams are not hard, perturbation theory is not reliable, and 
other descriptions of the corresponding reactions might be more adequate. If one treats them 
for instance as hadronic rescattering, there is again no reason why final state mesons X with 
small decay constants or higher spin should be suppressed. 

5.2 Annihilation contributions 

Annihilation diagrams are another important contribution violating naive factorization. In 
Tables ^| and |3| we have indicated the channels where they can occur, either from tree level W 
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(a) (b) (c) 

Figure 5: The different topologies for annihilation contributions to decays B — > YX, where X 
is one of our candidate mesons given in Table |I[ 

exchange or from penguin operators. An indication of the importance of annihilation could 
be obtained from data by comparing decays with and without such contributions that are 
otherwise similar, for instance B° — > D + a^ and B s — > D+clq , or B~ — > ti°D~j and B~ — > t]D~j. 

Depending on the flavor structure of the decay there are three annihilation topologies. The 
one shown in Fig. |5]a is for example relevant for B° — ► D + a$ and for B~ — > n^K^ . In 
this case only one of the quarks forming our candidate meson X originates from the effective 
weak vertex, so that our suppression mechanisms are not relevant. Notice that this holds true 
irrespective of whether the interactions between the gg-pair and the quarks attached to the 
decay vertex are under perturbative control or not, a point that is controversial in the literature 

HI- 

The two remaining topologies correspond to Zweig forbidden contributions. Neglecting 
electromagnetic interactions, they require that the meson formed from the gg-pair has isospin 
1 = 0. In Fig. |5|b, relevant for our decays into charmonium, neither of the quarks from 
the annihilation vertex enters in X, so that our suppression mechanisms again do not apply. 
Finally, there is the case of Fig. [5]c, where X is formed from the quarks of the decay vertex. 
Our suppression mechanisms only apply here if the gg-pair interacts solely with the quarks in 
the B meson, but not with the ones forming the meson X. 

The topologies for pseudo-annihilation contributions in b baryon decays have already been 
shown in Fig. |^. All of them circumvent suppression. 

6 The case of QCD factorization 

In the QCD factorization approach developed by Beneke et al. || [7], those corrections to 
naive factorization that are dominated by hard gluon exchange are calculated in perturbation 
theory, all other contributions are found to be power suppressed in 1/my,. The physical mech- 
anism underlying these results is color transparency of the meson ejected by the effective weak 
current. QCD factorization can therefore not be applied to decays where a D meson with 
its highly asymmetric quark- ant iquark configurations is emitted from the current, and we will 
thus not consider the corresponding channels in this section. 

As QCD factorization relies on color transparency, it also requires the emitted meson to be 
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fast in the B rest frame. Certain types of power corrections will therefore increase with the 
mass mx of the emitted meson. One can reduce the bias due to such effects by comparing our 
suppressed decays with unsuppressed channels involving mesons of similar mass. Examples are 
the p(770), ax (1260), or p(1450) for isospin-one mesons, and the K* and K\ resonances for the 
strange sector. 



6.1 Distribution amplitudes 

The most important effect of radiative corrections in our context is that the currents like 
(X\ d(0)7 M (l — 75)^(0) |0) occurring in naive factorization become nonlocal. The leading con- 
figurations in 1/rrif, involve light-like separations z and are parameterized by meson distribution 
amplitudes. In light-cone gauge, they read for a du meson 

(X(q,J 3 = 0)\d(z)Yh 5 H-z)\0) 2 = -ig" [' ^e^" 1 ) «>(«) + ... (25) 







to twist-two accuracy, where the . . . stand for terms of twist three and higher, which contribute 
to hard processes only at the power correction level. The Dirac matrices 7^ and 7^75 are 
to be taken for mesons with natural and unnatural parity, P = (— 1) J and P = (— 1) J+1 , 
respectively. The variable u gives the momentum fraction carried by the quark in the meson 
X, a natural frame of reference in our case being the B rest frame. Notice that the twist-two 
distribution amplitudes involving the vector and axial currents select the polarization state with 
zero angular momentum J3 along q in that frame. We remark that our subsequent discussion 
remains valid if instead of a meson X one considers a continuum final state with appropriate 



quantum numbers as discussed in Sect. In this case ip is to be replaced with a generalized 
distribution amplitude JEJ, defined as in Eq. ( ^B] ) with the appropriate replacement of state 
vectors \X). 

One easily sees that the lowest moment of ip in u gives back the local currents, so that for 
mesons with spin or 1 one recovers the decay constant, 

1 dui f (u) = f x . (26) 







The nonlocal currents in Eq. ( p5|) have however matrix elements for mesons of any spin J. 



Taylor expanding the bilocal operators around z = gives in fact operators with arbitrarily 
high numbers of partial derivatives between d(0) and u(0), and thus with arbitrarily high spin. 
In other words the lowest moment / du<p(u) of the distribution amplitude, projecting out the 
local V or A current, vanishes for mesons with J > 2, but not the function <p(u) itself. Hence 
the production of mesons with spin 2 and higher is no longer forbidden at the level of a s 
corrections to naive factorization. We thus find an explicit realization of our arguments in 
Sect. |5TT| : gluon exchange such as in Fig. f| indeed makes the production of higher-spin mesons 
possible. 



Let us now see what becomes of the other suppression mechanisms we discussed in Sect. pTT 
Charge conjugation invariance implies that the distribution amplitudes ipx for the neutral 
ao and bi mesons are odd under the exchange of quark and antiquark momenta, (px°(u) = 
—(px°0- — u), so that their first moment (|26|) is zero. In the exact isospin limit, the distribution 
amplitudes for the charged and neutral mesons in an isotriplet are the same, so that with 
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Eq. (EBT) the decay constants of the charged ao and b\ have to vanish, as already seen in 



Sect. |3J]. In the real world, the part of the distribution amplitude that is even under u — > 1 — u 
is therefore small for the charged ao and b\, i.e., ip x -(u) + (fx-(l — u ) ~ m d — iTiu- This does 
however not restrict the odd part of (fx- i u ), which can be comparable in size to the distribution 
amplitudes of, say, the 7r or the p. Considering the distribution amplitude of the Kq and using 
SU(3) symmetry we see that its even part is suppressed by m s — m u , but not its odd part. 

Along the same lines of reasoning, we find that up to isospin breaking effects the distribution 
amplitudes for the charged heavy pions are even under u — > 1 — u. The lowest moment 
/ du ^(1300) (u) is small of order m u + md according to Eqs. fll2|) and (p6p, but there is no such 
restriction on higher even moments such as / du (2u — l) 2 ^(1300) ( u ) ■ We thus see that neither 
spin nor any of our other suppression mechanisms applies at the level of a s corrections. 

It is useful to expand the distribution amplitude in Gegenbauer polynomials C^ 2 , which 



are the eigenf unctions of the leading-order evolution equation p3[ for quark distribution am 



plitudes. In order to achieve a uniform notation for mesons with different spins we write 

/i) = r 6u(l - u) [B + £ BM C 3 J 2 (2u - 1)] , (27) 

71=1 

where we have explicitly displayed the dependence on the factorization scale \x. For mesons X 
with spin or 1, we have Bq = 1 and f v is just the decay constant fx, whereas for mesons 
with J > 2 we have Bq = 0. We will in the following only use products f^Bn so that we need 
not specify the separate normalizations of f v and B n in that case. From our above discussion 
it follows that for our candidate mesons one or more of the lowest coefficients in the expansion 
( ^7) are either zero or small of order — m u , + m u , or m s — m u . In the Appendix we will 



estimate the orders of magnitude of the leading coefficients to be 

l/^iL^a^iq « 75MeV, 
\rB 2 \ <im ^ 2jP3 w 50MeV, (28) 

evaluated at the renormalization scale \i — m& = 4.4 GeV. 
6.2 Decays into a D and a light meson X 

For these channels the only diagrams one needs to consider at leading order in 1/mb are vertex 
corrections such as in Fig. 0a. The result of the 0(a s ) calculation for the matrix element of 
the effective weak Hamiltonian can be written as the sum of 

(D + X-\H cS \B°r rr = -i^V cb V: d a c rrq,(D + \crb\B°), 

(D* + X-\H eS \B°r™ = t^V cb V: d arrq^D* + \cYl 5 b\B°), (29) 

and the contribution of naive factorization, taken with a coefficient a^ act which equals the 
coefficient a\ of Sect. || evaluated at next-to leading order, up to a small term removing the 
renormalization scheme dependence. For details we refer to Eqs. (95) and (96) of 0. The 
coefficients a\ orr are given by 



rar = ^ CM ^F 1 duF ( ( (go) 

Ait l\ r Jo 
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a? 


(2) 


(3) 


(4) 
°1 


B° - 


-> D+X- 


12.3 + 12.7i 


-8.1 + 17.8i 


1.7-0.8i 


-2.3 + 4.9z 


0.3-0.5z 


B°- 


-> D*+X~ 


11.7 + 8.4? 


-9.5 + 16.0* 


0.3 - 1.2i 


-3.2 + 5.0i 


-0.2 -0.2i 



Table 5: The first few coefficients af 1 ' of the expansion fl3"T|) in units of 10 3 . They are 
evaluated at renormalization scale n = m b = 4.4 GeV with input parameters = 220 MeV 
and m c /m b = 0.3. With the same parameters one has a^ act = 1.039. 

where the function F(u,z) with z = m c /m b can be found in 0. Its second argument is +z 
for decays into D and — z for decays into D*, with z = m c /m b , so that af rr depends on both 
mesons in the final state. Its dependence on the distribution amplitude of X can be expressed 
as 

oo 

ar r (/i) = i?o + e 4 n V) £»(/*)■ (si) 

71=1 

The first few coefficients c4 are listed in Table |5]. We see that they are small compared with 
1, and that they tend to decrease with n. They depend substantially on the renormalization 
scale /i. One finds 

4 n) (m b /2) _ a ( r\m b ) _ ^ 
4 n) (m b ) 4 n) (2m b ) 

a falloff mostly due to the Wilson coefficient C%- The Gegenbauer coefficients B n also decrease 
with the factorization scale, although by less than a factor 1.2 for B\ and B2 when fi is varied 
between m b /2 and m& or between m b and 2m&. The effects of this dependence on /1 are quite 
mild for decays where most of the result is due to the Born level term /^a^ act , but not for our 
decays where this contribution is absent or suppressed by a small value of f v . A more stable 
prediction would require the inclusion of 0(a 2 s ) corrections. Since they involve again the color 
allowed Wilson coefficient C\ they may actually not be small compared with the 0(a s ) terms. 

Whereas a^ act is real valued, we observe from Eq. and Table [| that a™ rr is complex. The 
strong phases of unsuppressed channels like B° — > D + n~ are thus small in QCD factorization. 
On the contrary, they can be sizeable in our suppressed decays, where a s contributions are 
essential. 

6.3 Decays into charmonium 

The radiative corrections to factorization for decays into charmonium states have not been 
calculated yet. The following observation is however relevant in our context. The naive 
factorization formula for these decays involves the color suppressed coefficient a 2 and color 
suppressed penguins, but at the level of loop corrections the color allowed coefficient C\ will 
come in. Hence the 0(a s ) terms will probably be sizeable compared with the naive factorization 
result. This expectation is supported by an analysis of inclusive B decays into charmonium 
p3j] . Whereas naive factorization forbids the decays of Table |3] into a Xco or a Xc2 instead, one 
may then expect that within QCD factorization their branching ratios are not much smaller 
or even of similar size than for the corresponding decays into J/^f or Xci- 
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6.4 Decays into two light mesons 

For these channels, not only vertex corrections need to be considered, but also so-called penguin 
contractions (cf. Fig. 7 of 0) and hard interactions with the spectator quark from the B as 
shown in Fig. [|b. The latter involve the twist-two distribution amplitudes ipy and fx of both 
final state mesons. Mesons with J > 1 have a second twist-two distribution amplitude p5 



in- 



volving the nonlocal tensor current d(z) a^u(-z), which describes states with helicity J3 = ±1 
and can also contribute now. 

Beyond the level of leading twist contributions, several power corrections have been consid- 
ered in the literature || 26]. A particular type of correction that is numerically not suppressed 



occurs for final-state mesons with spin J = and involves their twist-three distribution ampli- 
tudes defined from the nonlocal (pseudo) scalar and tensor currents [E7|]. Their contribution to 
the amplitude relative to the twist-two radiative corrections is controlled by the ratio r defined 
as in (p2"D, which is formally of order l/rut, but numerically not small. With the quark mass 
values we use, one has r n — Tk — 1 5 and the corresponding twist-three radiative corrections 
have size similar to the leading-twist ones. This is also the case for the ao, 7r(1300), and Kq 
mesons we are considering here. Indeed, the twist-two distribution amplitudes are comparable 
in size for ao, 7r(1300), Kq, and n, K, as follows from comparing our estimates (|28|) with f n 
or fx- The same holds for the respective twist-three distribution amplitudes, which are con- 
trolled by the local (pseudo) scalar matrix elements in Eqs. ( |19D and (p0|). We recall that the 
twist-three pieces just discussed can only be estimated in QCD factorization because they lead 
to logarithmic divergences at the endpoints of the distribution amplitudes. 

Logarithmic endpoint divergences also appear in annihilation contributions, which in the 
power counting scheme of QCD factorization are again 1 jm^ corrections, even for those terms 
involving only twist-two distribution amplitudes. In a recent study of decays B — > tttt and 
B — > 7rK Beneke et al. [B] estimated annihilation contributions to be moderate corrections 
to the leading terms computed in the QCD factorization framework, giving a benchmark 
number of 25% in the branching ratio, although with large uncertainties. Notice that the 
small parameter controlling the relative weight of annihilation and leading contributions in 
their calculation is 

JbJy 



{m 2 B -m 2 Y )F^ Y {m\y (33) 

which depends only mildly on the mass of the emitted meson through the B — > Y transition 
form factor. We take this as an indication that the importance of annihilation contributions 
is not primarily driven by the size of mx- 

We expect then that in decays such as B° — > clqCIq the hard nonfactorizable terms calculable 
in QCD factorization, and possibly also the power corrections estimated there should be small 
compared to the amplitude for B° — > 7r + 7r~, which is dominated by the contribution of the 
large coefficient a\. 

For penguin dominated decays like B — > ttK, the overall size of corrections found in Ref. 
is not so small compared with the result of the naive factorization formula. With the "designer" 
modes B° — > -n^K^ and B~ — > n^K^ we can isolate such nonfactorizable terms. Among 
these, annihilation contributions from scalar penguin operators are of special interest. In QCD 
factorization they are a power correction, but not in the PQCD approach, where Keum et 
al. found that they contribute at order one and with a large phase to the B — > nK decay 
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amplitudes. Comparing the branching ratios of the above decays into K 2 with those into a K 
or K* can show whether scalar penguin annihilation contributions are indeed large. 



7 Branching ratio estimates for decays B — > DX 

We present now our numerical estimates of the branching ratios for decays B° — > D + X~ 



where X is one of a , a 2 , &i, 7r(1300), n 2 , P3, or _K"q(1430), K 2 . As discussed in Sect. |6T2j , these 
channels receive hard gluon corrections to naive factorization, which can be calculated in the 
QCD factorization approach. 

We recall that for the a 2 , vr 2 , P3, and K^-, the tree term proportional to a^ act is absent. The 
contribution from the a s correction a c ° TT given in (^) involves a contraction of with the 
matrix element parameterized as 



(D(j/)\crfb\B(p)) = F.iq 2 ) { (p + pY - 




+ F ( 



m, 



m 



V, 



(34) 



where q = p — p'. One easily sees that the a s contributions always pick up the form factor F , 
independent of the spin of X. Thus, the decay rate for our candidates X, as well as for any 
other spin zero meson like the 7r, can be written as 



T(B Q 



D + X~ 
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m 



2 \2 

d) 



rn 



■px 



B 



P 



,fact 



Fo(m 2 X/ 



(35) 



where px denotes the magnitude of the three-momentum of X in the B rest frame. 

We normalize the rate of the decays into X to the unsuppressed one into a light pion. This 
has the advantage that the CKM factors cancel (except for the strange mesons) and that we 
can use the measured branching ratio for B° — > D + n~ decays, where naive factorization works 
well B, 0] . For the ratio of decay rates we have the simple expression 



T(B C 



D+X- 



T(5 C 



D+ 



7T 



la^ + a? orr ) 



„fact 

7T U 1 



(36) 



where for simplicity we have neglected the term a™ rr for the ir~ , where its effect is of order 
1% to 2%. Eq. ( P6"D has corrections due to phase space and the evaluation of the form factor F Q 
at different momentum transfer q 2 , which go in opposite directions. The relevant q 2 ranges from 
roughly 1 GeV 2 for the ao(980) to 3 GeV 2 for the 7T 2 and P3, to be compared with q 2 for the 
7r. We have checked that the relation ( |36| ) is affected by these mass effects by not more than 
10% to 15%. This is sufficient in our context, given the dominant theoretical uncertainty of 
our calculation hidden in the decay constants and distribution amplitudes. For X = ^(1430), 
K2 we have to include a CKM factor | V^. s /V^d| 2 , which is known to a good precision. In our 
numerical analysis we take |V^ s /\4d| = 0.23 from flTT|. 



We evaluate ([36|) using the expansion ( |3T| ) with the coefficients of Gegenbauer moments 
in Table |5], the estimates (^) for the Gegenbauer moments B n , and the maximal values of 
the decay constants in ( |T4D to (0). The b\ can have a small decay constant, for which we are 
not aware of any information in the literature, and in our calculation we set it to zero. Before 
presenting the branching ratios let us study the magnitude of the individual terms entering in 
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We have /Uf ct = 136 MeV for the pion, and for the mesons with leading Gegenbauer 
moments f^Bi and f ip B2, we respectively obtain 

Lm<»,ks,k 5 = 1-5 MeV, 
|/^i° rr U(i300),. 2 ,p 3 = 0.1 MeV, (37) 

at renormalization scale u = mi,. The second term is tiny mainly because of the small coefficient 

(2,) 

a\ from the one-loop calculation. Since our estimate in the Appendix does not yield the sign 
of f^Bn, we have a twofold ambiguity when adding f^af" to f^a^, and find 

/ a o(98o) af* = 2.6 MeV, \f*> (af 1 + af) U(980) = 2.4 or 3.5 MeV, 

/ao(i450)< ct = 0.7 MeV, |H< rt + OU(i450) = 1.3 or 1.9 MeV, 

/.(1300) «l act = 7.5 MeV, \f* (a^ + a?") | ff(1300) = 7.4 or 7.6 MeV, 

/^(i430)« f i act = 48 MeV, |H«i act + < rr )k *(i430) = 47 or 48 MeV. (38) 



Notice the enormous correction to the naive factorization result for the a (1450). The im- 
pact of nonfactorizing corrections is similarly strong for the ao(980) if we take the minimum 
value / ao (980) = 0.7 MeV from Eq. (|i~5|). For the tt(1300), nonfactorizing terms remain always 
moderate since |/ v 'a5 orr | 7 r(i3oo) is small even compared with the lowest estimate of ^(1300) in 
Eq. (PS|). The K£ has a decay constant much larger than \f lp a c i rr \K* and is the only case where 
the corrections hardly matter. 

Using B(B° — > D~7i + ) = (3.0±0.4)-10 -3 from 0] and the maximal values in (pSf), we obtain 
the branching ratios in Table [!| To show their dependence on the choice of renormalization 
scale, we give them for /1 = rrib and fi = rrib/2. We take the latter as an indication of how 
large the branching ratios can be in QCD factorization, although they are not upper bounds 
in a rigorous sense. We also show the corresponding results from naive factorization, where for 
consistency of comparison we have again neglected the corrections discussed below Eq. fl36|). 
Here the scale dependence is minute, less than 2 percent, and the values in Table ^| are those 
for fi = rrib- 

We proceed to decays into a vector meson, B° — > D* + X~. The contraction of g M with the 
matrix element {D*\c r )^ r )^b\B) depends again only on a single form factor, commonly referred 
to as A and defined e.g. in ||. The decay rate is then given by 

r(B°^D* + X-) = ^p\\V c ^:/\r(a[^ + ar)A (m 2 x )\\ (39) 



As mentioned in Sect. |6.2j , the coefficients a™ rr for decays into D* are different from those into 
D, and we now have 

\r(a[ act + aP) L(980) = 2.2 or 3.5 MeV, 
|r(«i act + « C r r )laO(i450) = 1.2 or 1.9 MeV, 
ir(< ct + a C i° rr )U(i300) = 7.5 MeV, 
|r(4 act + ar rr )k*(i430) = 47 or 49 MeV (40) 

at fi = rrib. The analogous expression for the ratio (|36| ) of decay rates still holds, although 



with different mass corrections. We estimate them to be not much larger than for the D and 
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■ 10~ 6 


2.0 
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■ io- 8 


5.8 


io- 7 
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io- 7 
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10~ 6 


B° - 
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IO" 6 
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-> L>*+a (980) 


1.0 


• 10~ 6 


1.8 


10~ 6 


3.7 


IO" 6 


B° - 


-> L>*+a (1450) 


7.9 


■ io- 8 


5.2 


io- 7 


1.9 


IO" 6 


B°- 


-> £>*+&i, D*+a 2 







2.9 


io- 7 


1.5 


10~ 6 


B° - 


-> L)* + tt(1300) 


8.3 


■ 10~ 6 


8.4 


10~ 6 


8.4 


10~ 6 


B°- 









5.7 


io- 10 


3.2 


io- 9 


B° - 


-> L>*+ff*(1430) 


1.8 


■ 10~ 5 


1.9 


IO" 5 


1.9 


IO" 5 


B° - 


-> D* + K* 







1.5 


io- 8 


7.7 


io- 8 



Table 6: Branching ratios for various decay modes obtained in QCD factorization for two 
choices of the renormalization scale p. For comparison we also give the branching ratios in the 
naive factorization approach. We recall that we have set f bl = for lack of better knowledge. 



neglect them as before. Normalizing to B(B° — > D*~n + ) = (2.76 ± 0.21) • 10~ 3 from we 
obtain the branching ratios in Table which are somewhat smaller than the corresponding 
ones for decays into the D. With the exception of decays into K%,iT2, or p 3 , all estimated 
branching ratios are larger than 10 -7 and within experimental reach at the 5-factories. We 
expect similar branching ratios for B s decays into and the same candidate mesons X, up 
to SU(3) breaking effects. 

To facilitate comparison of the branching ratios of Table | 
a K, we estimate the latter as 



into K, 



and K% with those into 



B(B C 



D + K~ 



f 2 V 2 
Ik v " 



J IT 



V, 



ud 



f 2 V 2 
B(B° -> D*+K-) ~ J -§ £ 



ud 



B(B° 
B(B° 



D + 7i- 
D* + tt~ 



which is in good agreement with data on the ratios B(B° 
0.079 ± 0.011 and B(B° D* + K-)/B(B° -> D* + n-) = 0.074 ± 0.016 reported in |28 



(2.4 ±0.3) ■ 10" 4 , 
(2.2 ±0.2) • 10" 4 , 
D + K-)/B{B° -> D+7T- 



(41) 



A considerable source of uncertainty in our decay rate estimates are the unknown meson 
decay constants and distribution amplitudes. Let us illustrate this for the channel B° — > 
D + a (980). With the minimal value of / ao (980) in Eq. we obtain \f v (a\ iee + af n )\ a0 ( 9 m) = 
1.3 or 1.9 MeV at p = m b . The larger of the two possibilities corresponds to a branching ratio 
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B(B° — > D + ao(980)) = 5.8 • 10~ 7 , about three times smaller than the corresponding one in 
Table |. 

We have already seen the sensitivity of our results on the choice of renormalization scale. 
This uncertainty is largest for the branching ratios which are zero in naive factorization, with 
a variation by a factor of 5 to 6 between /x = m;, and fi = How important it is for the 

other channels depends on the actual size of the decay constants and distribution amplitudes. 
We also remark that the reliability of the QCD factorization approach for our "light" mesons 
with masses in the range from 1 to 1.7 GeV might be questionable, at least unless finite- mass 
corrections can be taken into account. 

Most important, however, is that the hard nonfactorizing contributions are small on the 
scale of the amplitudes for unsuppressed decays. With all uncertainties discussed above we find 
that / </3 a™ rr is at most a few MeV, i.e., less than 5% of /^af 10 *. One may well expect that soft 
corrections (or annihilation graphs when they can occur) are bigger than the hard ones. For 
all our meson candidates except the Kq they would then lead to considerably larger branching 
ratios than we have estimated. 

The decays into Kq are different in this context. Here the perturbative corrections to 
the prediction of naive factorization are quite small and their uncertainties less relevant, and 
further soft corrections might or might not overshadow the factorizing piece. Once the decay 
constant of the Kq is known experimentally, one should of course refine our estimates by taking 
into account the phase space and form factor corrections in (|3^) . 

7.1 Background from B° decay 

We see in Fig. |2] that the same final state of our signal mode B° — > D + X~ just discussed can 
be produced in the decay of the CP conjugated parent meson, B° — > D + X~ . As mentioned 
in Sect. [|, this background is CKM suppressed with respect to the signal. On the other hand, 
the signal mode is punished by a small decay constant, whereas the background goes with 
/d ~ 200 MeV. One expects the background-to-signal ratio to be large, since at the amplitude 
level X 2 fo/ fx ~ 0(1). We recall that no such background exists in decays into a strange final 
state, B° -> D+X; or B s -> D+X~. 

The background to B° — ► D + X~ can of course be removed by flavor tagging. Experimental 
discrimination between B and B is however challenging in decays with branching ratios of less 
than 10~ 5 , and it is worthwhile to see how far one can go without a flavor tag. Let us therefore 
investigate in more detail the branching ratios of the background decays for the case of the ciq. 
We parameterize the matrix element for the B — > a® transition in terms of form factors F£ 
and Fi as 

(ao(*/)|u7^&|B(p)> = Ffitf) |(p + pT ~ — ^ + W) — ^ ( 42 ) 
Assuming naive factorization, the decay rates for B° — ► D^ + Qq decays can be written as 



r(i"->A1 = t£ v * , ao/ Poo VcdKb M D y\F^ D )\\ 



16n m\ 
G 



2 



2 



T(B°^D*+a ) = ^vl a VJf^\ ai f D *y\F?{ml»)\ 2 , (43) 



23 



where a± is the universal coefficient for color allowed decays in naive factorization, introduced 
in Sect. 0. Using a\ = 1.03 we find 



B(B° -> D+a (980)) = 2.1 ■ 10 



/ \V V* 1 x 2 

-6 I I ^cd^ufe 



7.3 -lO" 4 / \ 200 MeV/ V 0.5 / 1.55 ps 




B(B° - D*+a (980)) = 1.9 ■ 10 6 ( 11^51 V f /d * V( g(<0 y_^_ (44) 
V 0V ;; ^7.3-10- 4 / V230 MeV/ \ 0.5 / 1.55 ps' V ; 

where we have indicated the sensitivity to several input parameters which at present have 
significant uncertainties. In particular, very little is known about the form factors for B — > X 
transitions. Chernyak JL6] has recently estimated the form factor _p 1 (0) B ^ a °^ 1450 - ) ~ 0.46 with 



light-cone sum rules, indicating only a small enhancement over the corresponding one into a 
pion, where he cites Fi(0) B ^ 7T ~ 0.3. That the form factor for B — > should rather be larger 
than the one for B — > n is also plausible in the Bauer-Stech-Wirbel approach ||29|| . To see this, 
consider the constituent qq wave function of a charged a . In the case where the q and q spins 
couple to S 3 = 0, it has a zero at momentum fraction u = ~ due to charge conjugation, up 
to small isospin breaking effects. Being normalized to one, this wave function is then more 
pronounced towards the endpoints u = and u — 1, and thus can have a greater overlap with 
the asymmetric wave function of the B than the pion wave function can. Further information 
may be obtained in relativistic quark models pO . 



We wish to point out that experimental information on the form factor F±(m 2 D ) can be 
obtained from semileptonic decays at q 2 = m 2 D , 

dT(B° 4 £-j? e ) _ _G|_ 3 ly |2| F a/ „2X|2 /.gx 
dq 2 - 2Alx zPa \V ub \ \* x {q )\ , ^40j 

where we expect similar statistics as for the semileptonic decays B — > ir,p with branching 
ratios of a few 10 -4 , if the form factors have comparable size. One may then relate Fq with 
F" using large energy effective theory (LEET) originally introduced in Ref. p2[ . We are 



in the kinematical situation where a light meson (here the ao) is emitted from a heavy parent 
with large recoil q 2 = m 2 D m 2 and an energy E = (m 2 B — q 2 + m^)/(2mg) much larger 
than Aqcd and the light masses in the process. This is the region of applicability of LEET. 
To leading order in 1/E and l/m^, we derive 

Fx a (E,m b ) = ^F a (E,m b ). (46) 

This is the analog of Eqs. (104) and (105) in pi] , to which we refer for details. Hence, the 
form factors are equal to leading order in the large energy limit. 

From the branching ratios in Eq. (]44]) we conclude that rather likely the background from 
decays of B° mesons into a does not overshadow the signal. A similar discussion can be 
given for decays into the other 1 = 1 mesons of Table [l]. Assuming naive factorization and no 
anomalous behavior of the relevant form factors, we quite generally expect branching ratios 
of the background modes B° — > D^ + X~ of order 10~ 6 . Any significant excess over both 
this and the branching ratios given in Table ^| would imply that either there are important 
nonfactorizing contributions in the signal, or that naive factorization drastically fails in the 
background channel. 
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8 Summary 



We have explored how to obtain quantitative information on nonfactorizing effects in exclusive 
b decays, using channels where such contributions are not hidden behind larger factorizing 
pieces. We achieve this through "switching off" the factorizing contribution by choosing final- 
state mesons with either a small decay constant or spin J > 2. Our proposal is similar in 
spirit to the study of decay channels where the quark content of the final state does not admit 
factorizing contributions, such as B° — > K + K~ [[33|], B s — > tt + tt~, 7r°7r°, or 6-decays into 
baryon-antibaryon pairs (see e.g. @ and references therein). Suppression of the factorizable 
contributions thus highlights factorization breaking effects, such as annihilation graphs, soft or 
hard interactions, and in general any mechanism dominated by long-distance physics, which 
disconnects the b decay vertex from the final state meson. We have explicitly shown that hard 
nonfactorizing contributions, calculated in the QCD factorization framework, can yield sizeable 
contributions to the decay amplitude. 

In a systematic study, compiled in Tables ^|, ^|, and|], we have shown that our method applies 
to a variety of mesons and channels, in decays of B u ^, B s , and b baryons. In particular, the 
mesons X we have selected cover a wide range of masses, which makes it possible to explore 
whether the energy- mass ratio Ex/mx m the parent rest frame is a relevant parameter to 



ensure factorization, as is suggested by color transparency but not by large N c arguments p2 



We have presented a detailed analysis of color allowed decays B — > D^*'X and B s — > DyX 
with a light meson X. When the factorizable contribution is suppressed, e.g. for the scalar 
ao, we found that hard nonfactorizing corrections can be of similar magnitude or even larger 
than the Born term. They remain however much smaller than the amplitudes of corresponding 
nonsuppressed decays, for instance into a tt. In several cases we found branching ratios sub- 
stantially enhanced over the ones calculated in the naive factorization approach, see Table |6|, 
and are within the reach of existing and future experiments at the S-factories BaBar, Belle, 
CLEO and at hadron colliders like the Tevatron and the LHC Comparison of these decays with 
modes where the ejected meson is a D meson and further study of those into charmonia Xco or 
Xc2 should give complementary information on the origin and limitations of the factorization 
approach. 

B decays with light-light final states are more complex. Modes such as B° — > a^a^ and 
B s — > K + ciq are not entirely suppressed by the small decay constant of the a due to the 
presence of scalar penguin operators, but we find the corresponding amplitudes to be much 
smaller than those of B° — > 7r + 7r~ or B s — > K + n~. Such factorizing penguin contributions 
can be eliminated altogether with higher-spin mesons like the b\ or 02 instead of the ao- We 
expect hard nonfactorizing contributions to be moderate, too, so that experimental information 
on such decays could again tell us whether nonperturbative effects are large. For penguin 
dominated decays like B° — > -k + K^ the situation is less clear-cut on the quantitative level, 
but we argue that they can give valuable indications on the importance of penguin annihilation 
contributions. This is a particularly controversial issue since different conclusions regarding 
such decays have been drawn in the QCD factorization and the PQCD scenarios || ||. 

To conclude, we find that the b decays presented here provide a tool for studying impor- 
tant issues in exclusive nonleptonic decays. We stress that in order to make this tool more 
quantitative, the decay constants of the ao(980), a (1450), 7r(1300), _K"q(1430) and b\ mesons 
should be known experimentally. Their determination from r decays should be in reach of the 
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existing experiments at BaBar, Belle and CLEO, and even more of dedicated r-charm facto- 
ries. Information on the distribution amplitudes of ao(980), ao(1450), 02, 7r(1300), 7C2, which 
are needed for the calculation of hard nonfactorizable contributions, could be obtained from 
7*7 collisions at the B factories. 

Note added: The suppression mechanisms discussed in this work provide opportunties to 
study CP violation in various B decays into designer mesons a , &i, 0,2, etc. This has been 



explored independently in two recent studies 37 and 38 . 
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Appendix 

In this appendix we estimate the size of the leading-twist distribution amplitude p of several 
mesons. Our method is based on the connection between distribution amplitudes and the Fock 
state expansion in QCD, and closely follows the discussion in |[34j| , to which we refer for details. 



The starting point is to decompose a hadron state on Fock states consisting of current quarks 
and gluons, qq, qqg, etc. The coefficients in this expansion are the light-cone wave functions 
for each parton configuration. For a du meson, one has 

\x ) J3=0 - / i u[l _ u) m^ 71 V>(«,*l)+..., (47) 



where u and k± denote the light-cone momentum fraction and transverse momentum of the 
d quark in the meson. The arrows indicate quark and antiquark helicities, and the + and 
— respectively apply to mesons with natural and unnatural parity. The states \d^ui) and 
\di Mf) are understood to be coupled to color singlets. By ... we have denoted the Fock states 
\d^u{) and with aligned quark helicities, and Fock states with additional partons. The 

connection of the light-cone wave function ip(u, k±) with the distribution amplitude defined in 
Eq. (H) is 

,hkl ^( u ,k ± ) = -^= V (u). (48) 



167T 3 rV ' ^' 2V6 

The probability to find the du Fock state with antialigned helicities in the meson X is 

du^\i:(u,k L )\ 2 . (49) 

This should be below 1 since it is the probability to find a current qq pair in the meson, without 
further gluons or sea quark pairs. Note that this is different from the qq wave functions in 
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constituent quark models, which are by definition normalized to 1. Let us now use the relation 
( f49|) to estimate the size of <p(u). In order to achieve this, we need to make an ansatz for the 
k± dependence. A form consistent with several theoretical requirements |34], |35 



is 



ip(u, k ± ) = -j^- a exp 
u(l — u) 



a 2 k 2 



u(l — u) 



2^6 



(50) 



where the prefactor of the exponential is imposed by the normalization (|48|). a plays the role of 
a transverse size parameter of the dw-pair in the meson. For the pion, Brodsky and Lepage 1531] 
obtained « 0.86 GeV -1 with the above ansatz and the asymptotic form <p n (u) = f w 6u(l— u) 
for the distribution amplitude. This corresponds to an average transverse momentum (k\) ~ 
(370 MeV) 2 and to a Fock state probability of P^ ~ 0.25. We will take the same values for 
the mesons we discuss here, which is certainly a crude assumption but should give the correct 
order of magnitude. 

With the ansatz (|50f) and the Gegenbauer expansion (|27| ) we obtain 



P = 2n 2 (ap 



, g 3(n + 2)(n + l) g2 



n=l 



2(2n + 3) 



(51) 



Consider now a meson for which Bq = 0, such as the ai or K\- If we take a = a T and 
P = P n ~ 0.25 and retain only the term with £?i in the Gegenbauer expansion, we obtain 



\FB 1 



100 MeV. 



(52) 



Including the zeroth term J^Bq in the Gegenbauer expansion, as is appropriate for the charged 
ao, Kq and 6i, would decrease this estimate by about 5% for the Kq when taking fx* = 42 MeV. 
The effect of that term for the ao or b\ can be neglected even more safely. 

In order to explore the dependence of our estimate on the ansatz we made for the k± 
dependence of the wave function, we take an alternative form 



ip(u, k A 



16vr 2 5 4 



u 



k\ exp 



~a 2 k 2 



u 



2VQ 



<p(u), 



(53) 



which has a node at k± = 0. For the Fock state probability we find the same expression as 
fl5T|) with a replaced by a/\/2. Choosing a = a we then get an estimate of l/^-Bi) larger by a 
factor y/2. If instead one requires the average (k'jj to be the same with the two forms ( pop and 
(|53D, one finds a = \/3a and thus an estimate of [f^Bx] smaller by a factor of y^2/3. Given 
these observations we expect that fl52"D should give the correct order of magnitude of the first 
Gegenbauer coefficient. 

We should add that this does not hold for mesons that are not qq bound states in the 
constituent quark picture but for instance made from qqqq, which may be the case for one 
of the ao mesons. Is plausible that for such a system the probability of finding a single 
current qq pair in this meson is reduced compared with the one of a conventional qq state. 
Correspondingly, its twist-two distribution amplitude <p and the coefficients f v B n would be 
smaller than estimated here. 
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Our considerations are easily adapted to the case of mesons where B\ is zero or isospin 
suppressed, such as the 7r(1300), the 7r 2 , or the ps. Retaining only B 2 in the Gegenbauer 
expansion of <f(u) and taking as before a ~ a n and P = P n m 0.25, we obtain 

\f v B 2 \ « 80 MeV. (54) 

So far we have not displayed the dependence of both the distribution amplitude <p and 
the light-cone wave function ip on the factorization scale fi, which physically represents the 
resolution scale of the qq pair. Our above estimates are understood as corresponding to a 
hadronic scale, say, fi = 1 GeV. Evolving up to \i — m b we obtain \f v Bi\ w 75 MeV from 
Eq. (ID and \f tp B 2 \ « 50 MeV from Eq. (|§. 

To conclude this section, we wish to point out that experimental constraints on the distri- 
bution amplitudes for the neutral mesons X = ao, a 2 , 7r(1300), tt 2 can be obtained from the 
process 7*7 — > X at virtualities Q 2 of the photon much larger than the meson mass. This can 
be measured in e + e~ — > e + e~X, and the CLEO data for X = 7r, 77, 7/ are in fact one of our best 
sources of information on the corresponding distribution amplitudes ||36|| . To leading order in 
l/Q 2 and in a s , the amplitude for 7*7 — > X is proportional to f v (B + B 2 + B± + . . .) for 
mesons with unnatural parity, and to /^(-Bi + B 3 + B 5 + . . .) for mesons with natural parity. 
According to our above estimates, one then expects cross sections comparable to the one for 
7r production, so that the measurement of these reactions at large Q 2 may well be in the reach 
of the B factories. 
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